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I. INTRODUCTION 


Various ways were tried recently to decrease the friction drag of 
a body in a flow; they all employ influencing the boundary layer 
(reference l). One of them consists in keeping the boundary layer 
laminar by suction; promising tests have been carried out by Holstein 
(references 2 and 3 ) and Ackeret (reference 4). Since for large Reynolds 
numbers the friction drag of the laminar boundary layer is much lower 
than that of the turbulent boundary layer, a considerable saving in 
drag results from keeping the boundary layer laminar, even with the 
blower power required for suction taken into account. The boundary layer 
is kept laminar by suction in two ways: first, by reduction of the 

thickness of the boundary layer and second, by the fact that the 
suction changes the form of the velocity distribution so that it becomes 
more stable, in a manner similar to the change by a pressure drop 
(reference 7) • Thereby the critical Reynolds number of the boundary 
layer (U5*/v) cr it becomes considerably higher than for the case without 
suction. This latter circumstance takes full effect only if continuous 
suction is applied which one might visualize realized through a porous 
wall. Thus the suction quantities required for keeping the boundary 
layer laminar become so small that the suction must be regarded as a 
very promising auxiliary means for drag reduction. 

Various partial solutions exist at present concerning the theoretical 
investigation of this problem. Thus H. Schllchting (references 5 and 6 ) 
investigated the plane plate in longitudinal flow with homogeneous 
suction. At large distance from the leading edge of the plate a constant 
boundary layer thickness and an asymptotic suction profile result. Later 
H. Schlichting and K- Bussmann (reference 9) investigated the two- 
dimensional stagnation point flow with homogeneous suction and the plate 

in longitudinal flow with v Q ~ 1 /fx (x = distance from the leading edge 
of the plate) . In all cases a strong dependence upon the mass coefficient 
of the suction resulted for the velocity distribution and the other 
boundary layer quantities. K. Bussmann, H. Munz (reference 8 ), and 
A Ulrich (reference 16) calculated the transition fran laminar to turbulent 
(stability) of the boundary layer with suction for several cases; in all 
of them the stability limit was found to have been raised considerably by 
the suction. As is known from earlier investigations (reference 7), the 
same amount of influence on the transition from laminar to turbulent is 
exerted by the pressure gradient along the contour in the flow for 
impermeable wall. Both influences (pressure gradient and suction) will 
be present simultaneously for the intended maintenance of a laminar 
boundary layer for a wing. Both influences have a stabilizing effect for 
the suction in the region of pressure drop; in the region of pressure rise, 
however, presmire gradient and suction have opposite influences. Whereas 
without suction, for pressure rise, mostly transition in the boundary 
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layer occurs; here the important problem arises whether this transition can 
he suppressed by moderate suction. 

The solutions for the laminar boundary layer with suction existing 
so far are not sufficient for answering these questions. An exact 
calculation of the boundary layer with suction encounters insuperable 
numerical difficulties .Just .as in the case of the impermeable wall. Thus 
it is the more important to have an approximation method at disposal 
which permits one to check the calculation of the boundary layer with 
suction for an arbitrary body. Such a method will be developed in the 
present treatise. The method given here is an analogon to the well-known 
Pohlhausen method for impermeable wall. It permits the calculation of 
the laminar boundary layer with suction for an arbitrarily prescribed 
shape of the body and an arbitrarily prescribed distribution of the 
suction velocity along the contour In the flow. 


x,y 


5* 





II. SYMBOLS 
(a) Lengths 


coordinates parallel and perpendicular, rec actively , to 
the wall wetted by the flow (x * y * 0: stagnation 

point and leading edge of the plate, respectively) 


displacement thickness of the boundary layer 

momentum thickness of the boundary layer 

N 

u/u(l - u/br)a y 

S 

measure of the boundary layer thickness 
plate length or wing chord, respectively 
plate width 



(1 - u/u) dy 1 



(b) Velocities 

u,v velocity components in the friction layer, parallel and 

perpendicular to the wall 

U(x) potential velocity outside of the friction layer 

U 0 free stream velocity 

v 0 (x) given suction velocity at the wall v 0 < 0: suction 
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Q 
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(c) Other Quantities 


wall shearing stress 


wall shearing stress for asymptotic solution of houndary 
layer on plate in longitudinal flow with homogeneous 
suction 


total suction quantity 



dimensionless mass coefficient of suction; Cn > 0: suction 

/uSn 


reduced mass coefficient of suction 


(cq 



q dimensionless distance from wall (y/Sq) 

F^(rj), ^(tj) basic functions for velocity distribution in boundary layer, 

equations (8), (9) 

K form parameter of boundary layer profiles, equation (6) 

dimensionless boundary layer thickness, equations (12), (13) 
dimensionless momentum thickness, equations (22), (23) 

I dimensionless length of boundary layer 

III. THE EQUATIONS OF THE BOUNDARY LAYER 

WITH SUCTION 



Following we shall consider the plane problem, thus the boundary 
layer on a cylindrical body in a flow (fig. l) . x, y are assumed to be 
the coordinates along the wall and perpendicular to the wall, respectively, 
U Q the free stream velocity, U(x) the potential flow outside of the 
friction layer, and u(x, y) , v(x, y) the velocity distribution 
in the friction layer. Suction and blowing is introduced into the 
calculation by having along the wall a normal velocity v Q (x) prescribed 
which is different from zero and generally variable with x: 


v Q (x) > 0: blowingj v 0 (x) < 0: suction 

v o/ u o ma y assumed- to be very small ( 0.01 to 0 . 0001 ). Only the case of 
continuous suction will be considered, where, therefore, vo(x) is a 
continuous function of x. One may visualize this case as realized by a porous 
wall. The tangential velocity at the wall should, for every case, equal 
zero. The boundary layer differential equations with boundary conditions 
are for the steady flow case 
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B 2 u 

By^ 


( 1 ) 


Bu Bv 
Bx By 


0 


(la) 


y = 0: u 

y =<» : u 


0; v 
U(x) 



( 2 ) 






The system of equations (l), (2) differs from the ordinary 'boundary 
layer theory merely "by the fact that one of the boundary conditions for 
y = 0 is changed from v = 0 to v = v Q (x) ^ 0. Thereby the character 
of the solutions changes decisively: the solutions differ greatly 

according to whether it is a case of v Q > 0 (blowing) or v 0 < 0 (suction). 


A special solution of the equations (l), (2) which forms the basis 
for the theory of the boundary layer with suction and is used again 
below is the solution for the plane plate in longitudinal flow with 
homogeneous suction, thus v Q (x) = v Q = const < 0 and U(x) = U 0 . Tor 
this case the boundary layer thickness becomes constant at scone distance 
from the leading edge of the plate} also, the velocity distribution 

becomes independent of x (reference 5) • From = 0 follows because 

Bv 

a? 


of the continuity 


r 0 and hence 


v(x, y) = v Q = Constant 


From equation (l) then follows for the velocity distribution 
u(x, y) = u(y) = U 0 


( ^r) r -y/W\ 

- e J = U 0 ^l + e J 


( 3 ) 


with 


= 


-V r 


(M 



6 


NACA TM No. 1216 


signifying the displacement thickness of the asymptotic solution. The 
wall shearing stress for this solution is: 

T o = s “P u o T o (^a) 

W/ 0 

It is independent of the viscosity. This asymptotic solution is one of 
the very rare cases where the boundary layer differential equations can 
be integrated in closed form. 

For solution of the boundary layer differential equations (l), (2) 
for the general case where the contour of the body and hence U(x) and 
also v 0 (x) are prescribed arbitrarily one could consider developing 
the velocity distribution from the stagnation point into a series in 
terms of x in the same way as for impermeable wall (reference 11); 
the coefficients of this series then are functions dependent on y for 
which ordinary differential equations result. K. Bussmann (reference 17) 
applied this method for the circular cylinder; with a very considerable 
expenditure of time for calculations the aim was attained there. However, 
for slenderer body shapes the difficulties of convergence increase so 
much that this method which works directly with the differential equations 
is useless for practical purposes. 


IV. THE GENERAL APPROXIMATION METHOD FOR 
ARBITRARY PRESSURE DISTRIBUTION 
AND ARBITRARY DISTRIBUTION OF 
THE SUCTION VELOCITY 
(a) The Expression for the Velocity Distribution 


For that reason one applies an approximation method which uses 
instead of the differential equations the momentum theorem which represents 
an integral of these differential equations. By integration of the 
equations (1), (2) over y between the limits y = 0 and y = « one 
obtains in the known manner (reference 18 ): 


U 2 — + (26 + 6*) U— - Uv = T — 
dx dx o o 


(5>) 


j 
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■6 signifies the momentum thickness, 6* the displacement thickness, and 


calculation of the boundary layer to he chosen here proceeds in such a 
manner that a plausible expression is given for the velocity , distribution 
in the boundary layer u(x, y) which is contained in equation (5) 
in S* and t 0 . Thus an ordinary differential equation for ^(x) 

results from equation (5)» after this differential equation haB been 
solved one obtains the remaining characteristics of the boundary layer 
&*(x) , T 0 (x), and the velocity distribution u(x, y) in the boundary 

layer- The usefulness of this approximation method depends to a great 
extent on whether one succeeds in finding for u(x, y) an expression 
by appropriate functions. 

Pohlhausen (reference 15) first carried out this method for the 
boundary layer with impermeable wall. The velocity profiles in the 
boundary layer were approximated by a one -parametric family and the approxi- 
mation function for the velocity distribution expressed as a polynome of 
the fourth degree. The coefficients of this polynome are determined by 
fulfilling for the velocity profiles a few boundary conditions which 
result from the differential equations of the boundary layer. This 
method proved to be satisfactory for the boundary layer without suction. 

Thus one proceeds in the same way for the boundary layer with 
suction. For the velocity distribution in the boundary layer one chooses 
the che -parametric expression 



the wall shearing stress- The approximation method for 


U ■ *i(t>) + KF2(tj) > q 


y 

5 x (x) 


( 6 ) 


F-^(q) and Fg(q) are fixed prescribed functions which are immediately 
expressed explicitly; K = K(x) is a form parameter of the boundary 
layer profiles, the distribution of which along the length is dependent on 
the body shape and the suction law; Si(x) is a measure for the local 
boundary layer thickness. The connection between 8q and 6* and is 
given later. It proved useful to choose other expressions for the 
functions Fj_(q) and F 2 (q) than Pohlhausen for the impermeable wall. 

For the velocity profile according to equation (6) the following five 
boundary conditions are prescribed; they all follow from the differential 
equations of the boundary layer with suction, equation (l), (2): 

(7a, b) 


y.O: u-0; v 0 |=u| + v|| 


sy 

6u 


o 2 u 


u = U Q ; ^ = 0; ~ = 0 


^y 


£y 2 


yaco: 


(7c,d,e) 
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The selection of (tj ) and. F 2 (q) is to "be made from the view point 

that a few typical special cases of velocity profiles of the "boundary 
layer with suction are represented "by equation (6) as satisfactorily as 
possible- In particular we shall require the asymptotic suction profile 
according to equation (3) to be contained in the expression (6). This 
condition is satisfied if one puts 


=» 1 - e _T 1 


( 8 ) 


and correlates the values K = 0 and Sj_ = 5* to the asymptotic suction 
profile. Furthermore, the expression (6) naturally should yield usable 
results also for the limiting case of disappearing suction. To this 
purpose a good presentation of a typical boundary layer profile without 
suction Is required. One chooses as this profile the plate flow for 
impermeable wall according to Blasius (reference 11) • Since no convenient 
analytical formula exists for the exact solution of this case, a good 
approximation formula for Blasius' plate profile is needed. It Is found 

that the function = sin(ari) gives a very good approximation to the 
Blasius profile (a = Constant )A Thus one puts 


0 < ti < 3: F 2 (b) * Fi - sin(|i) 
■n 3 : F2 (t) = F-j_ - 1 = -e _T l 



and then obtains with K = -1 a good approximation for the plate flow 
without suction. The corresponding value of &]_ is given later. The 
functions Fi(t|) and F 2 (t]) are given in figure 2 and table 1. Thus 
one has for the velocity distribution in the boundary layer the expression: 


x That the sine function is a good approximation for the velocity 
distribution at the plane plate without suction, resulted from an 
investigation of Mr. Iglisch about the asymptotic behavior of the plane 
stagnation point flow for large blowing quantity (reference 2o) . 
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0 < n < 3 : 


li = 1 - e _T W K |jL - e _T > 



q > 3: ^ = 1 - (K + l)e"H J 


( 10 ) 


By selection of the functions Fi and F2 the "boundary conditions (7a,c,d,e) 
are per se satisfied. The last boundary condition, equation (7b), results, 
because of 


_o 

P 


*s V 


si [ 1 + 4 - £>] 


(13) 


in the following qualifying equation for K: 

If [1 ♦ <1 - 1 )] - UU' -vJL(l.K) 


and from, it with 


Bq^U' 

X - - 7 — 

h" 

ro 

, " v o 5 l 

X 1 * V 

(13) 


for K the equation 


X + Xq - 1 
1 - - €> 


(U0 


\ and Xq are two dimensionless boundary layer parameters. A quantity 
analogous to X was already used by Pohlhausen for the boundary layer 
without suction*, Xq is newly added by the suction. For the asymptotic 
suction profile with 8 q = 5, Xq = 1 according to equation (4). The 
form, par ame ter K as a function of X and Xq is represented in 
figure 3 . 
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(b) The Differential Equation for the Momentum Thickness 

In order to obtain by means of the expressions (6), (8; and (9) 
from equation ( 5 ) the differential equation for the momentum thickness 
one must first set up the relations between 3,5* and S-,. For the 
displacement thickness there results: 



(1 - Fi)drj - K 



(15) 


The calculation of the integrals gives: 


- = 1 - K<2 - D * g*(K) 


(16) 


For the momentum thickness one obtains: 


r“ f 
61 Jo 


(Fi + KF 2 )(1 - Fi - KF 2 )cLti 


§£ = C Q + C X K + C 2 K 2 = g(K) 


(IT) 


The calculation of the integrals gives: 


“ J F]_(l - Fi)dr| = ^ 


(18a) 


Ci - 


* 00 

f (F 2 - -F^dr, = 


n 


-1. f - * — 

« 3 ,^2 


6e3 


1 + 


(D‘ 


0 


= 0.06656 (18b) 
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t 


C 2 = “I F 2 2 dt| 


- 3 + ~ 
J n 


1 + 

n 6e^ 
3 

1 + 


= -0.02358 (18c) 


Hence 


C 1 - 0 2 = 2 - -2 = °. 090U 


Thus there Is 


ST - | + 0.0665QC - 0.02358c 2 = g(K) 


(17a) 


For the form parameter of the 'boundary layer profiles &*/8 used later 
one obtains therefore 


8^ 

d 


1 - g(g - 1) 

^ + C^K + C2K 2 


1 - Q.Q9014K 
+ 0.0665QC - 0.0235& 2 


(19) 


Furthermore there results according to equations (n) and (IT) : 


^ S’ - g [l + k(i " I)] = f(K) (20) 


The functions g(K) , 6*/S and T o&i/dU according to equation (11) are 

represented in figure h and table 2. 

In order to derive from equation (5) a differential equation 
for d(x) one writes equation (5) in the form 


Ud dd 



jo * 
u u 


( 21 ) 


Furthermore one introduces according to Holstein and Bohlen (reference 19) 



12 


NACA TM No. 1216 



and 


-v Q -d 

— “ K 1 = X lS 


( 22 ) 


(23) 


With 


Z = 


*2 


v 


then 


k = ZU'j 



(24) 


(25) 


is valid. With equations (22) to (25) as well as equation (20) the 
differential equation (2l) is transformed into 




+ icj_ » f(K) 


(26) 


If one finally puts for abbreviation: 


G(k , K^) = 2f - 2k 


1 - K 


2 + 


C 2 -D 


g(K) 


- 2K, 


(27) 


the differential equation for Z(x) becomes: 


dZ G( k, ki) 
dx “ U 


K 



( 28 ) 
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If the function G(k, *ci) is known, the integral curve Z(x) 
can "be calculated from this equation by means of the isocline method. 

For carrying out the calculation in practice it is useful to 
introduce dimensionless quantities. One forms them with the aid of 
the free stream velocity U 0 and a length of reference l (for instance, 
chord of the wing) . Thus one puts 


Z* 


Zll 


_o 
l > 


x * - 



- f x (x*) 


(29) 


Then equation (28) becomes: 


d _ G(ie, *l) , 
dx* U/U 0 * 



dU 

dx> 


K 1 = f“i(x*) /z* 


(30) 


The function g(k, «i) is calculated as follows: First, one obtains k 

and as functions of X and Xx from equations (22) and (23), if 

one takes the connection between K and X, X^ according to equation (lM 
into consideration: 




K B g 2 (K)X = g 2 (x, X;l)X 
K j_ = g(E)Xj_ = g(X, Xi)Xi 

From equations (27) and (31) follows: 


(31) 




1 2 

1 - =(? - D 

I - \(T 

2 + 

_ s 


iis 


|g = g^ 


! + k(i - - 2Xg - x[l - - D] - X-H 


G - 2gF(X, X x ) 


(32) 
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with 


F(X, X x ) = 1 + K(l - -g)- 2Xg - X [l - K<2 - |)] - 


(33) 


Hence G can "be calculated first as function of X, Xj_ and then, 
■because of equation (3l)> also as function of k, k-j.* 

The functions k(X, Xj.) and Kj_(X, X^) are represented in figure 5 
and table 3* The function thus determined G(k, k^) is given in 
figure 6 and table 3* 


(c) Stagnation Point and Separation Point 

The behavior of the differential equation (28) at the stagnation 
point where U = 0 requires special considerations. In order that the 
initial inclination of the integral curve (dz/dx) Q _ at this point be 
of finite value, G( k, must equal zero. This g’ives the corresponding 

initial values Kq, k^ q . Since the function g(K) does not have a zero 
for the values of K considered (compare fig. 4) the determination 
of the initial values k 0 * amounts to the zeros of 

F(X 0 , X 1o ) * 0 (34) 


The resulting initial values at the stagnation point X Q , Xj 0 are given 
in table 4, together with the initial values Kq, calculated 

additionally according to equation (3l). To each pair of values K o> B lo 
corresponds a mass coefficient of suction which results from 



-v 0 (o)d 0 

*o» v = *lo 


as 

-Vq(o) *10 _ 

/ “ y — c °0 

F'o v / B o 

In figure 7 the initial values and *2 0 80:6 plotted against the 

local mass coefficient at the stagnation point. The initial value Z c 
corresponding to k q is obtained by 
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K 


Z„ » 


° ' Uo' 


(35) 


The following connection exists "between the distribution function of 
the suction fi(x) = and C Q : Uq' * KiU 0 /i, being 


a profile constant. Thus there is 


and 



f l(°) 

ft 


(36) 


The determination of the initial values of the integral curve proceeds, 
therefore, as follows: With the given initial value of the suction 

velocity at the stagnation point f^o) one first determines C Q 
according to equation ( 36 ). One obtains the corresponding initial 
values * 0 and K ^o from figure 7, and according to equation ( 35 ) the 
initial value Z Q of the integral curve. If the suction does not 
begin at the stagnation point but further downstream, 0 o » 0 ; 

*lo - 0 and according to figure 7 


K 0 = 0.0709J Z c 


0.0709 

V 


Separation point . - The separation point is defined by the fact that 
the wall shearing stress there equals zero. This gives for K, according 

to equation (ll), the value K = = -2.0 99* For the asymptotic suction 

* 6 - it 

profile K = 0; this is simultaneously the greatest possible value of K. 

To K = -2.099 corresponds the value \ = —- 1 - = -1.099 for all and 

6 - n 
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k = -0.0T21 for all *!• However, if one would want to carry the 
■boundary layer calculation up to this point, certain difficulties result 
in the last part shortly ahead of this point, since the correlation 
between it and \ is not unequivocal there (compare fig. 4) . The 
function g(k, ki) against k also is not unequivocal shortly ahead 
of this point. Thus it is useful to select a point situated somewhat 
further upstream as separation point where the boundary layer calculation 
has to stop. Such a point results if one chooses the K - value of an 
exact separation profile according to Hartree (reference 13)* For this 
latter there is: 


separation: % 



- 0.0682 


(37) 


One. defines this point as separation point of the present boundary 
layer calculation for all mass coefficients of suction. The following 
table gives a survey of the values of -8 and S* at the separation 
point for four different calculation methods: 


Case 

fv' > K k 

, . * 
v u = X A 

5* 

3 

New method: (sine 

approximation ) 

(-0.0721) 

(-1.55) 

(4.64) 

Pohlhausen P4 
(reference 15) 

-0.1567 

-1.92 

3*50 

Exact Hartree 
(reference 13 ) 

-0.0682 

-1.22 

4.03 

Exact Howarth 
(reference 14) 

-0.0841 

-1.25 

3-84 


The selection of the separation point thus made is somewhat arbitrary; 
however, it may be accepted unhesitatingly since, as is well known, the 
approximation methods for the boundary layer calculation in the region 
of the pressure rise are always somewhat uncertain and only a rough 
estimate but no exact calculation of the boundary layer parameters is 
possible here- For the same reason one may also accept the fact that 
for the present case the velocity distribution u/U 0 partly assumes, 
shortly ahead of the separation point, values which are. slightly 
larger than 1. 
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(&) Performance of the Calculation for the General Case 


By means of the system of formulas given above one may perform the 
calculation of the boundary layer for an arbitrarily prescribed body 
shape and an arbitrary distribution of the suction velocity along the 
wall in the flow. It takes the following course: 


To the distribution of the suction velocity -v Q (x) corresponds 
the total suction quantity 


Q = 


v Q (x) dx = “CQU 0 bi 


x=0 


and the reduced mass coefficient 


CQ* = 



( 38 ) 


and the reduced suction distribution function according to equation (29) : 


fl(x*) 


~Vq(**> ,/Pq* 

u 0 I V 


Thus there is 


c q * = 


f x (x*) dx* 


c *=0 


( 39 ) 


If the suction begins at the stagnation point, one determines 

with Kq = ^ Qbc) the mass coefficient C Q at the stagnation point 

according to equation (36). Then one obtains k q and ^ rcm 

figure 7 and Z Q according to equation (35)* With these initial values 
the differential equation (30) can now be graphically integrated by means 
of the diagram in figure 6. The calculation is carried out up to the 
point where k reaches the value = "0.0682. This integration 
immediately yields Z*, k, Kq as function of x*, with 
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* 

if v 


(bO) 


The remaining "boundary layer parameters then result as follows: 

By means of figure 5 one obtains after k and the parameters X. 
and snd additionally from figure 3 the form parameter K. After K 

one obtains from figure 4 the form parameter 6 */t 3 and thus 



From equation (20) one then also obtains the wall shearing stress t q 


Jo fn f(|0 

nuf u 0 = fz * 


(40a) 


Finally, the parameter is required for the velocity distribution in 

the boundary layer. According to equations (6), (7) > and (40) : 




I ,/UflJ g(K) 

* If v & 


(^1) 


Examples of such boundary layer calculations are given in chapter VI. 


V. SPECIAL CASES 
A. WITHOUT SUCTION 


Our general system of formulas is to be specialized in this section 
for a few typical special cases for which one can partly give solutions 
in closed form. First, the case without suction in particular shall be 
treated for which, naturally, our equations also must give satisfactory 
results. This case one obtains for v D (x) = Oj then 
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\ ± 5 0; K]_ s 0 (without suction) 

and equation (14) is transformed into 

K = X - 1 (42) 

Therewith, according to equation (17): 

g^- = g(x) “ 5 + C^(x - l) + C 2 (l - l) 2 

■ + 0.06656 (x - 1) - 0.02358 (X - 1) 2 (43) 

The differential equation (28) for the momentum thickness "becomes 


= SLSl; K = zu' (44) 

dm u 

G(k) Is, according to equation (32) and (33) : 

0 = 2gF(x) (45) 

F(X) = 1 + (X - l)(l - |) - 2X |j + C X (X - 1) + C 2 (X - l) 2 ] 

- X [l + (1 - X)(2 - D] 

F(X) = -2C 2 X 3 + (2C 2 ' 2 + f) X 2 (l " | - |) *■ + (46) 


Furthermore, according to equation (3l) : 
it = g 2 X = X 


r + C^(X - l) + C 2 (x ~ l) 


'] 


(47) 
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The values of G and k calculated according to equations (43), (45), 
(46), (47) are given in table 3- 


(a) The Plane Plate in Longitudinal Flow 

The boundary layer at the plate in longitudinal flow without suction 
(Blasiue) wich U = Uo is obtained for X = k = 0. Then, according 
to equations (42) to (46): 




K = -1 


g(0, 0) * \ - c x + C 2 = 


i . 1 
« 2 




F(0, 0) = - 


G(0, 0) * 2|(| - |) = 2 - | = 0. 


429 


J 


(48) 


With the initial value Z Q = 0 the integration of equation (44) then 
gives: Z-(5-f)i or 



For the form parameter 8*/$ 


follows from equation (19): 



2.66 


4 - n 


(50) 
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and thus for the displacement thickness 


S* = (it - 2) 


k - « r U, 


¥ - 1 - 740 ,/^ 


U, 


( 51 ) 


For the coefficient of the total friction drag Cf 


Vi 


R 


iu 0 2 hl 


of the 


plate of the width h and the length l, wetted on one side, one 

28 i 

obtains because of Cf. = — — 

v 

1-308^ (52) 

Finally, the velocity distribution is, according to equation (10) , 


0 < ri < 3 : u = U 0 sin 

ti > 3 : u = U Q 


Therein Is ^ 



and 





(%) 


In figure 8 the velocity distribution according to equations (53) 
and (51) is compared with the exact solution of Blasius} the agreement is 
very good. Furthermore the characteristics of the boundary layer according 
to the present approximation calculation are compared with the values 
of the exact solution of Blasius in the following table. For further 
comparison the values according to the approximation method of Pohlhausen 
(reference 15) also have been given- The agreement of our new approxi- 
mation method with the exact solution is excellent for all boundary layer 
parameters} the drag coefficient, in particular, shows an error of 
only 2 percent. 
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Coefficients of the Boundary Layer at the Flat Plate 
in Longitudinal Flow Without Suction 


Calculation method 

let: 

* 

to 

,/s 

f vx 

&* 

d 


t 0 d 

n u 0 

New method (sine 
appr oximat i on) 

1.740 

0.655 

2.66 

1.310 

0-215 

Pohlhausen P4 
(reference 15) 

1.750 

0.685 

2.55 

1-370 

0.234 

Exact (Blasius) 

1.721 

0.664 

2-59 

1.328 

0.220 


The deviations of our sine approximation from the exact solution are, 
for most characteristics, even somewhat smaller than in the Pohlhausen 
method. 


(b) The Plane Stagnation Point Flow 

For the plane stagnation point flow the velocity of the potential 
flow U(x) = u-jx. All boundary layer characteristics are in this case 
Independent of the length x. The initial value of the momentum 
thickness Z 0 is obtained from equation (44) for G(k 0 ) = 0. Since g(K) 

does not vanish in the range of the values of K considered, there 
must be F(X, 0 ) = 0. From equation (46) one finds as zero of F(\) the 
value 


\ a = 0.3547 (stagnation point without suction) (5*0 


The corresponding values of K and K according to equations (42) 
and (47) are K Q = -0.6453 end « 0 = 0.0709> furthermore there is, 
according to equation (43): g(\ 0 ) = 0.447. Therewith the momentum 

thickness for the plane stagnation point flow becomes: 


d 



1/ ,TT 


ui 


0.266 



(55) 


The form parameter 8*/d results from equation (19) as 5*/d = 2.37; 
therewith one has 
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&* - 2.3T/k7 0.630|/^ (56) 

Furthermore, according to equation (43) : 8 ^ = 0-595 Thus there 

results from equation (ll) for the vail shearing stress: 

S /I * 1 - 1 ® (57) 

The velocity distribution results from eq.uation (10) as: 


“\ 

0^. ti < 3 : JJ- - 0-3574(1 " e -T 0 + 0-6453 

TI 3 : ,f *=1 ’ 0-3574e _T l 

u o 

J 


(58) 


with 


*1 



( 58 a) 


Figure 8 gives a comparison between velocity distribution according to 
equation ( 58 ) and the exact solution by Hiemenz (reference 12 ) j here also 
the agreement is satisfactory- Furthermore the characteristics of the 
boundary layer according to the present calculation are again compared 
with the exact solution by Hiemenz and with the approximate calculation 
by Pohlhausen in the following table- 
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Coefficients of the Boundary Layer of the Plane 
Stagnation Point Flow without Suction 


Calculation method 

*/? 


5* 

■d 

is? rr 

nu l*i 

t 0 i3 

pu 

New method (sine 
appr oximati on) 

0.266 

0. 630 

2.37 

1.163 

0.310 

Pohlhausen P4 
(reference 15) 

0.278 

0. 64l 

2.31 

1.19 

0.331 

Exact (Hiemenz) 

O.292 

0.648 

2.21 

1.234 

O.36O 


The agreement of the new method with the exact solution is for this 
case somewhat less satisfactory than for the plane plate; neither is it 
quite as good as the approximation of Pohlhausen. But even here the new 
method yields still very useful values. 


B. WITH SUCTION 


In this section a few cases with suction will "be treated far which 
the solutions can he given In closed form. First we shall treat the 
boundary layer at the plate In longitudinal flow with homogeneous suction, 
already investigated formerly (reference 6) . The following results are 
considerably more accurate than those former ones. 


(a) Growth of the Boundary Layer for the Plate 
in Longitudinal Flow with Homogeneous Suction 

For this case the boundary layer is at large distance from the 
leading edge of the plate independent of x; hence all boundary layer 
parameters are constant. The corresponding asymptotic solution has 
been given already in equations (3), (4), (4a). The applying values 
are: 


& *co ■ ZFi F ■ 2 f T o„ ■ -pu 0 v 0 ; K = 0 


> 


6 1 ■ 6 1o » ^1*1} K 1 = 1 


J 


( 59 ) 
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One now calculates the growth of the "boundary layer from the value 
zero at the leading edge of the plate to the given asymptotic value. 

In our system of formulas one has to put for it: 


X — 0 j k — 0 


Therewith becomes according to equation (l4) with the abbreviation 


1 - ^ = c = 0.4764 


K = 


- 1 


1 “ c\]_ 


and according to equation (17) 


g- » g(0, X^ = 


Po + Pl^l + P2^l' 
(1 ' cXi) 2 


with 

p o “ ~\ + % = 0.40986 

p l " ^ + 1 ' 7 " I 3 ! - °' 06656 

*2 = - 2 + I + 2 ^ " if + ?Pl " °'° 5819 


(60) 

(a) 


(62) 


(62a) 


Furthermore there is according to equation (16) 


6 * 


3 - f + *1^3 + n + |) 


(63) 
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The vail 


shearing stress becomes 


according to equations (ll) and (13): 


T 


O 


Jt 

Uo 6 



b l 1 - 0\l 


-pv 0 U 0 


it 

6 

Xi(l - cXi) 


(64) 


The differential equation (28) assumes for the present case the form: 


dZ _ G(0, *l) . 
dx = U Q ’ 


*1 


-v 


o 



v Q =* Constant < 0 


(65) 


The integration of this differential equation requires the explicit 
expression for G(o, k^). According to (32) and (33) : 


G(o, ki) 


2g 


f-i- 

\1 - cX x 



< CA '1 

* 2g(l - ^2) 

1 " cXt 


( 66 ) 


Thus G(o, k^) ■ 0 for *» 1. Therefore X^ = 1 is a solution of 
the momentum equation; it corresponds to the asymptotic solution. The 
initial value at the leading edge of the plate is Xi = 0. For the 
length of growing boundary layer Xj varies from 0 to 1. 

If one Introduces as dimensionless distance along the plate 



(67) 


the differential equation (65) can be written In the form: 


d(*i 2 ) 

dl 


G{*l) 


( 68 .) 


Initial value: 1=0; = 0 


(68a) 
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The connection between tc x and. X x is given by 

% 

*1 “ *-is(\l) 


(69) 


with g(X X ) according to equation (62)* 


The differential equation (68) can be solved according to the 
isocline method. For the present case, however, an analytical solution, 
too, is possible which is preferable. From, equation (69) first follows: 


2*! 


dic x dl x 
dXq dS 


G(k x ) 


(70) 


Here all quantities can be best expressed by X x so that a differential 
equation for X x (£ ) results. With dic x /dX 1 according to equation (69) 
one obtains from equation (70) after division by 2g since the latter 
does not disappear in the range 0 ^ X x £ 1: 


i 



dg ''NdXi 

+ dx^aT 


(1 - x x ) 



(71) 


Initial value: 5=0: X x = 0 


(71a) 


Because of g(o, o) = p Q = ~ one obtains from it in the neighbourhood 

of 1 = 0, X 1 = 0: 


and 



(72) 


I = 


2 

Tt 



0.391X3 2 


(7?) 
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Hence follows for the neighborhood of the leading edge of the plate 
(£ = 0) 

Xi «= 1. 60 /i or * 1. 60 

v 


Because of 



follows hence: 


~ y o & * 

v 


( * - 2 ) 



or 

5 * ■ (« - 2> 

As the comparison with equation (51) shows, the boundary layer thickness 
starts, therefore, at the leading edge of the plate with the value for 
the plate without suction. 

In order to integrate the equation (71), one has to insert the 
explicit values of g(X.i) and dg/dXi according to equation (62). 

After same intermediate calculating (compare appendix i) one obtains 


with 


dl_ „ X lfco + P1M. + P2M 2 + P^l 3 ) 

6X1 (1 “ cl 1 ) 2 (l - ll)(_| - cx^) 


*0 

Pi 

?2 


Po = 0 . 40986 

2pi + Po c = -O.i+667 

3P2 = 0.17457 i 


P3 = -cp2 = - 0.02772 J 


(74) 


( 75 ) 
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The breaking up into partial fractions yields: 


41 ^ Ki Kg 

- + + + 




dXi c 3 - 1 c x x - | (cXi - l) 2 cX-2. “ 1 


f'Ul) (76) 


The integration -with the initial value X-j_ * 0 for 1=0 yields 


I = %X X + K L In <1 - Xj.) + In (l - fs X ^) 


Ko + ~ In (1 - oXi) = f(X X ) 

J clj. "1 c 


(77) 


The K-p • • . , K^ result from the breaking up into partial fractions as 


K l = -6.9560; Kg » 3.4704; K 3 » -0.2284; K k = -0.1569 (78) 


Thus the solution finally reads 


1 = -0.2564\ 1 - 6-956 In (1 - Xi) + 7-2846 In (1 - 0.9099XJ.) 

X»-i 

+ 0.228k — ! - 0.3293 In (1 - 0.4764x x ) (79) 

2 

For development of this solution in_ the neighborhood of I = 0, 

X 1 = 0, the coefficient of Xj_ must, because of equation (73), equal zero; 

the coefficient of X^ must equal ^ = O.39I. The result is: 

^ - Kl - fcg + K 3 - K k = 0 (77a) 


*1 

2 


lfo , r 

jT K 2 


+ cK- 


" I k 4 


2f6 . 

it 


(77b) 


With the numerical values of equation (78) one may verify that these 
equations are satisfied. 
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The solution calculated accordingly is given in table 5* 

From X]_(£) all remaining boundary layer parameters can then be calculated 


immediately according to equations (6l), (62), ( 63 ), and (64). They also are 

-v 0 &* 

given in table 5* In figure 9 — y — t 


8* 4 
t and 


T 0 8* 


HU 0 


are plotted against j/T. 


The displacement thiclmess of the boundary layer reaches 0. 95 of 
its asymptotic value after an extent of the growing boundary layer of 
2 

4.5* The velocity profiles in the growing boundary 
, -voy 

layer in the plotting u/U 0 against — ; — = are represented in 

figure 10* For the wall shearing stress one obtains from equations (64) 
and (4a) : 


Is_ 

t o m 


n 

z 


Xj(l - cX^) 


(80) 


The wall shearing stress is plotted in figure 11 as a function 
of /T- 

Drag. - In view of the reduction of the drag by maintenance of a 
laminar boundary layer the friction drag in the extent of growing 
boundary layer is of particular interest for this solution. For the 
asymptotic solution the local friction drag along the wall is constant 
with t O0# » -pv 0 U 0 j thus the coefficient of the total friction drag 
also equals this value 


c f. 


W 


T 0 


Ju o £ bl 


|» 0 


00 

2 




(a.) 


For small suction quantities -v 0 /u 0 the extent of growing boundary 
layer is sometimes so large that the growth is not finished by far at the 
end of the plate. According to former investigations (references 8 and 10) 
it is to be expected that for homogeneous suction at the plate the 
maintenance of a laminar boundary layer is possible even for Reynolds 

Uq v *~T A 

numbers of the order of magnitude nj- =10* to 10 with a very small 

suction quantity of the order of magnitude cq = = 10'^. For — = 10 

U 0 Z 7ft U ° U ° 

= 10 ' or 10° one has at the end of the plate 

/'v 0 \2 u n z 

( Uo ) v = 0.1 or 1, that is, the growth of the boundary layer 


and 

s 


is not finished by far. 
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Since the friction layer over the extent of growing boundary layer 
is thinner, the friction drag there is considerably larger than for the 
asymptotic solution. For this reason the calculation of the drag over 
the extent of growing boundary layer will be given completely. 

The total friction drag for the plate wetted on one side is: 


W = b / t q dx 


( 82 ) 


and with the value of T 0 according to equation (8o) and with T 0ec 
according to equation (4a): 


w j 0 ,Tirr^ 


2 

With dx ■ tf“ according to equation (67) this equation becomes: 


2 v * / d£ 

p ° 


(83) 


with signifying the value of 5 at the end of the plate, thus: 


l ’l "(S) “ f(X lo) 


(84) 


Therein Vj_ Q signifies the value of at the end of the plate which 

is obtained from, equation (77) for I = therefore 


fUlo) = “3 ^-lo + In (l - Xio) + £ ^ ^ ^1^ 


- K- 


TLo 


°4io “ 1 


+ ^ In (l - cXi 0 ) 


(85) 
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Introducing In equation (83) for = f'(Xi) the expression according 


to equation (76) one obtains: 


W - PW 0 2 “-F(Xio) 
* o 


and 


W 

c f B p 0“ 

|u o 2 b2 


2-^F(X l0 ) 

“V- b> 


(86) 


with F(Xio) signifying 


F(Xiq) 


it 

I 


() x 1o 


x l“° 


f'(xi) dXi 
X-j^l - 0X1) 


(87) 


Finally introducing — — r * ,-r- 2 f(Xi G ) according to equation ( 8^4-) into 

v 0* u o 

equation (86) one obtains 


■T n r(\io) -t 


* ■ V ® ■ V 0(Xl ° ) 


or 


c f ■ c f#e O(Xio) 


( 88 ) 


Because of the connection between X;j. 0 and Ij according to 
equation (8b), the total drag coefficient for the extent of growing 
boundary layer is thereby given as a function of the dimensionless 

.2 


distance along the plate 


( 1 1 - 


On the other hand, equations (8b) 
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% 


a pd (86) give for prescribed mass coefficient of the suction -v 0 /[ J Q 
the drag law cf also against U 0 l/ V in the form of a parameter 
representation. The parameter Xio is the dime nsionless boundary 


layer thickness at the end of the plate: 



The values 


of X lo lie between 0 and 1, the first value being valid at the 
leading edge of the plate, the latter for the asymptotic solution, after 
the growth of the boundary layer has ended* 


The calculation of the integral F(Xi 0 ) according to equation ( 87 ) 
gives (appendix II) 


* 


fF ( xio > ■= < 3 § - rri - T - k 3 * K t ) to(1 - 

-t- ln(l - X io ) + ^K 2 In (l - ^rX lc J) 


* 


- ( -K3 + K4) 


cXio k 3 c Xio( 2 “ cX -lo^ 
1 - cXio 2 (1 - cXio) 


(89) 


c is, according to equation ( 60 ), c » 1 - g*, and Ki, • ••, K 4 are 
given by equations (75) and ( 78 ). After insertion of the numerical values 
according to equations ( 60 ) and ( 78 ) follows: 


F(Xio) 58 “O .3288 In (l ~ O.k76kx,2o) “ 6,956 In (l _ Xio) 


+ 7,28h6 In (1 - 0.9099Xio) 
0.4764 x 1o 


- 0.037^ 


1 - O.k764x lo 


0.05980 


O . h 764 x lo (2 - q. 476^X 1o ) 
(1 - o.VT64x 1o ) 2 


The values of F(Xio) and G(Xi 0 ) are given in table 6. For Xi Q — >1 
that is ^ — >« (growth of boundary layer ended) one has, as can 
immediately be seen from equations ( 89 ) and (77): 
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Xj_q — V ” I * G(Xio) = 1 (89a) 

f (^lo) 

and thus Cf — 5 Cf^ for i j— On the other hand one has in the 

neighborhood of the leading edge of the plate, that is, for Xi 0 — > 0 
according to equation (73) 


Xio “>0: f(Xi 0 ) = l (| - £) X l0 2 


and thus according to equation ( 87 ) 


>•10 -*0i F(X lo ) - (| - £) Ho 


and therefore 


X lo — »0: G(X lo )=f~ (89b) 

j ^lo 


If one substitutes this value into equation (88) and takes into consideration 
that 


x lo 



7Wd 


is valid for small Xq 0 according to equation (73)? one obtains 
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or 


c f = (Y)' 1/2 


thus the drag law of the plate without suction according to equation (52)- 
The drag law of the length of growing boundary layer is therefore for 
very small lengths of growing boundary ij asymptotically transformed 
into the drag law of the plate without suction. 

The drag law according to equation (88) is represented in figure 12, 
where Cf/cf^ is plotted against £j. Furthermore figure 13 gives the 

drag law in the form Cf against U 0 l/v for various values of the 
mass coefficient -v 0 /u 0 . The larger the suction quantity the smaller 
the Reynolds number at which the respective Cf - curve separates from 
the drag curve of the plate without suction and is transformed, after a 

-2v q 

certain transition region, into the asymptotic curve Cj^ - ~ — • The 

Reynolds number at which the latter is reached is the larger, the 
smaller the suction quantity. 

The drag coefficients given here represent the total drag of the 
plate with suction. No special sink drag is added (compare reference 10) 
since for continuous suction, as in the present case, the sucked particles 
of fluid have already given up their entire x- momentum in the boundary 
layer so that this momentum is contained in the friction drag. 

In order to obtain the total drag power of the plate with suction, 
however, one must, aside from the drag given here, take into account the 
blower power of the suction. 


(b) The Plate Stagnation Point Flow with Homogeneous Suction 

Another special case which can be solved in closed form is the plane 
stagnation point flow with homogeneous suction. Since for this case the 
exact solution from the differential equations of the boundary layer has 
been given elsewhere (reference 9) it shall also briefly be treated here- 
The potential flow is U(x) = uqx and the suction velocity 
v 0 (x) = v 0 (o) = v OQ <0. If the integral curve of equation (28) is to 
have a finite value at the stagnation point x = 0, there has to be 
G( k, Kq) = 0; this in turn requires, as was discussed in detail in 
chapter IV c, F(X, Xq) = 0. F(X, Xq) is given by equation (33). The 
values of X Q , Xq Q which belong together follow from it; they give for 

the general case the initial values of the boundary layer calculation at 
the stagnation point; for the present case of stagnation point flow they 
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immediately give the complete solution since the "boundary layer thickness 
and all other parameters are independent of the length of growing boundary 
layer x. Besides Xq, X one further obtains K according to 
equation (lk), g according to equation (17)# &*/$ according to 

equation (19)# and k 0 and Kq 0 according to equation (31)* The mass 

coefficient C Q = 7=^== is obtained according to equation (3k). From. * 0 

finally follows d and therewith 5*. The results are compiled in 
table 4. Natur ally # momentum and displacement thickness decrease with 
increasing suction quantity. With C 0 — > » the form parameter 5*/d 

approaches the value 2 of the asymptotic suction profile. In figure lk 

d fijv and S* y/uq/v are plotted against C Q and compared with the 
exact solution. The agreement is quite satisfactory. 

As conclusion of these considerations of the special cases the 
characteristic boundary layer parameters for these special cases are 
compiled in the following table. 



r 


Boundary Layer Parameters for Various Special Cases 


CaBe 

K 

X 

X 1 

It 

*1 

5* 

6l 

to 

II 

1 

Plate without suction 

-1 

0 

0 

0 

0 

1.089 

0.410 

2. 66 

Stagnation point flow 
without suction 

-0. 6453 

0.3547 

0 

0.0709 

0 

1.058 

0.447 

2.37 

(Separation point 
calculated) 

-2.099 

-1-099 
for all Xp 

f(cQ) 

-0.0721 

f(cQ) 

1.189 

0.256 

4.64 

Separation point 

according to Hartree 

f(X, Xi) 

f(Xp) 

f(c Q ) 

-0.0682 

f(CQ) 

f (cq) 

f(°Q) 

4.03 

Asymptotic suction 
profile 

0 

0 

1 

0 

1 

2 

1 

A 

2 

2 


VI. EXAMPLES 3 


In this section the now method shall be tried on a few more examples. 

(a) Circular Cylinder 

As first example the circular cylinder with homogeneous suction has been calculated for 
various suction quantities C Q . The results (displacement thickness 6* and form parameter it) 

3 The numerical calculations of thlB section have been carried out by Mr. A. Ulrich. 


U) 
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are given in figure 15* For the case without suction cp = 101.7° results 
as separation point j this is slightly further to the front than for the 
customary Pohlhausen method. (<p = 108.9°) for which the calculation was 
performed elsewhere (reference 7) • With increasing suction quantity 
results a reduction of the boundary layer thickness and a shifting of the 
separation poirit toward the rear. 

In order to completely avoid the separation for the circular cylinder, 
it is probably useful to select not a homogeneous suction along the 
contour, as in the present case, but a distribution of v Q (x) which has 
considerably larger values on the rear than on the foreside. Euch 
calculations may also be carried out according to the present method 
without additional expenditure of time. 

A comparison of the present approximate calculation with an exact 
calculation by K* Bussmann (reference 17) for the displacement and 
momentum thickness is given in figure 16* The latter calculation is a 
development in power series starting from the stagnation point, as first 
indicated by Blasius (reference ll) • Except for the neighborhood of 
the separation point the agreement is quite satisfactory. 


(b) Symmetrical Joukowsky Profiles for c a = 0 

As second example a symmetrical Joukowsky profile of 15 percent 
thickness has been calculated far c a = 0, also with homogeneous 
suction. The suction extends over the entire contour. The same profile 
without suction has been calculated elsewhere (reference 7) , also 
according to the Pohlhausen method. Here, too, a reduction of the 
boundary layer thickness and a shifting of the separation point toward 
the rear results with increasing suction quantity. For the suction 
quantity C G = 0-417, that is fj_(o) = 3> a separation does no longer 
occur. 


71 I. SUMMARY 


A method of approximation for calculation of the laminar boundary 
layer with suction for arbitrary body contour and arbitrary distribution 
of the suction quantity along the contour of the body in the flow Is 
developed. The method is related to the well-known Pohlhausen method for 
calculation of the laminar boundary layer without suction. The calculation 
requires the integration of a differential equation of the first order 
according to the Isocline method. The method is applied to several special 
cases for which there also exist, in part, exact solutions: Plate In longi- 

tudinal flow and plane stagnation point flow with homogeneous suction. 
Furthermore the circular cylinder and symmetrical Joukowsky profile with 
homogeneous suction were calculated as examples. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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VIII. APPENDIXES 
APPENDIX I 

Concerning the Length of Growing Boundary Layer for the Plane Plate 

with Homogeneous Suction 


According to equation (71) is 



(1 ~ Vl) 



(71) 


From equation (6?) one finds: 



(1 ~ clj)(p 1 + gpgl-jj + 2c (1 - cXi)(p 0 + pi^i + pgXx 2 ) 

(1 - cXx)^ 


Pi + g cp 0 + (PxC + 2gg)X 1 
(l - cXi)' 


( 1 , 1 ) 


• Substitution of e<iuatians(l,l) and (62) into equation (71) gives: 



fr 

O 




Pr> + 


v j- 


Pl x l + P2^1 

" „;_\s — 




PX + 2cp Q + (pxc + 2pa)>-l 


\L' 




_ 


(1 - M.) 


6 ' cX l IS 


1 a A\ 1 

ih WV JL 


and after multiplication ty (l - cXx): 

CPo + Pl*l * PS^l g )( 1 ~ °^l) + *-l(Pl + g cPo) * 4i £ (j?i c + 2 P&) 

(1 - c^x) 2 


dS 


u-M)(f-=^ (1,2) 


(1 


T ^Jg [po + ( 2 Pl + c Po)^l + 3^x 2 " c P2 x 1 3 J = (1 - ll)(| - c\x) £“ 


d g _ *-l(*o + ^1 + W* + *&lp 
^1 (1 - cXx) 2 (l - \j) Q - cXj) 


with 


D H * V. 

*0 - *01 






P2 » 3 P 2 > P3 ■ "cpg 


i 


( 74 ) 


( 75 ) 


«- 


t 
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APPEHDIX II 

Concerning the Calculation of the Drag of the Plane Plate 
with Homogeneous Suction 


The calculation of 


F(Xlo) * 7 


* (l X lo f'Ui) 


X^O 


Xj_(l - cl£) 


gives with f'(X.i) according to equation (76), if is replaced "by z, 


6u,/, > ^ lo dz „ lo dz 

? 10 = ° 3 L - -> + 1 L.o 


+ K2 


=0 ( cz " g) 2 ^ 1 " cz ^ 3 (i z =0 (° z " l) 2 z(l - cz) 


" K k 7 775" = 1+11+111+17 + 7 

4 (j z _o 


( 11 , 1 ) 


The integrals are solved by breaking up into partial fractions. One 
finds: 


I = ■— lln z - In 


Q -D] 


-K, iln z + In (z - l) + In (z - — ) 

1 c - 1 1-c V QS 

*- -1 z=C 



k-2. 
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III = -K 2 


[l ln 2 + c 2n (* " c) " 4 In Q ‘ 5g)] 

*“ z=0 


IV = F 2n z - Zn(cz - l) - ' 1 ~ 


L 


^ (cz - l) 2 


Mo 


-* 0 


V = Ki 


r t i 

| -In z + 

1 cz - 1 


+ 2n(cz - 1) 


" *lo 

i 

- 0 


When summed up, all terms with 


[ !n i 1 


cancel each other, because 


of equation (77a). After insertion of the limits the remaining terms give: 


|ln(cz - l)J = 2n(l - cXi Q ) j Jln(z - l)j = ln(l - X^ Q ) 


Mo -cX ^ 0 


[» • • & r ■ ■»<; - [it+iF-: 


cX 


lo 


1 ~]Mo cX, lo^ 2 ■ cX, lo) 


l Mo 
(cz - 1)^_ Q 


(1 - cX lo )‘ 
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Thus there results fcy simplification from equation (ll,l) 


jffclo) = ^3 " 1^1 ' 1? " e 3 + *k) ln(l ■ cX lo) 


+ 


Mi - Mo) 



cl lo Ko c?t lo^ 2 “ c ^lo) 
1 “ CX^q (l “ c ^>lo)^ 


” (-Kg + Kij.) 
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APPENDIX III 


To page 30 - " Table 10 gives the numerical table concerning the 
velocity distribution of figure 10. 

To page 3^ . " For the boundary on the plate in longitudinal flow 
with homogeneous suction the exact solution from, the differential 
equations also was given in an unpublished report by Iglisch. A 
comparison of the approximate solution above with that exact solution 
is given in figures 18 and 19* Figure 18 gives the comparison of the 
displacement and momentum thickness; particularly for the displacement 
thickness the agreement is good. Figure 19 gives the comparison for the 
wall shearing stress; here also the agreement is satisfactory. (So far, 
these comparisons can be carried out only for the front part of the 
length of growing boundary layer, up to I = 0. 5, since the exact 
solution does not yet completely exist.) 
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TABLE 1 


TEE BASIC FUNCTIONS AND F2 

FOE THE VELOCITY DISTRIBUTION 
IN THE BOUNDARY LAYER 
WITH SUCTION 


s) 

II 

F 1 

*2 

0 

0 

0 

.2 

.1813 

.0768 

.4 

.3297 

.1221 

.6 

.4512 

■ 1423 

.8 

-5507 

.1452 

1.0 

.6321 

.1322 

1.2 

.6988 

.1112 

14 

• 753^ 

-.0843 

1.6 

.7981 

•0551 

1.8 

.8347 

.0263 

2.0 

.8647 

-.0017 

2.2 

.8892 

-.0240 

2.4 

• 9093 

-.0416 

2.6 

.9257 

-.0524 

2.8 

• 9392 

-.0552 

3*0 

.9502 

-.0498 

3*5 

♦9698 

-.0302 

4.0 

•9817 

-.OI83 

4-5 

.9889 

-.0111 

5-0 

■9933 

-.0067 

6-o 

• 9975 

-.0025 

7-0 

•9991 

-.0009 

00 

1 

0 
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TABLE 2 


PARAMETER OE BOUNDARY LAYER WITH SUCTION 


-K 

e-i- 
8 &i 

* 5 * 
6* - 8J 




a 0 

0.5 

1 

2 

1 

0.5 

.1 

.4931 

1.009 

2.05 

.9524 

.4696 

.2 

.4856 

1.018 

2.10 


.4393 

•3 

.4779 

1.027 

2.15 

.8571 

.4096 

.4 

.4$96 

1.036 

2.21 

.8094 

.3801 

•5 

.4608 

1.045 

2.27 

.76L8 

• 3510 

* * 6 

.4516 

1.054 

2.33 

• 7142 

.3225 

13 .6453 

.4472 

1.058 

2.37 

.6926 

• 3097 

>7 

.4419 

a. 063 

2.41 

.66® 

.2945 

.8 

.4317 

1.072 

2.48 

.6L89 

.2672 

mm 

.4210 

1.081 

2.57 

• 5712 

.2405 

□jttM 

.4099 

1.090 

2.66 

.5236 

.2146 


.3982 

1.099 

2.7 6 

.4760 

• 1895 


.3862 

1.108 

2.87 

.4283 

.1654 

■EH 

.3736 

1.117 

2.99 

.3807 

.146l 

1.4 

.3606 

1.126 

3-12 

• 3330 

.1201 

1-5 

.3471 

1.135 

3-27 

.2854 

• 0991 

1.6 

•3331 

1.144 

3.43 

.2378 

.0792 

1-7 

.3186 

1-153 

3-62 

.1901 

.0606 

1.8 

.3037 

1.162 

3-82 

.1425 

.0433 

1-9 

.2883 

1.171 

4.09 

.0948 

.0273 

2.0 

.2724 

1.180 

4-33 

.0472 

.0129 

<*2.099 

.2562 

1.189 

4. 64 

0 



0 


Asymptotic suction profile. 
"Stagnation point without suction. 
c Plane plate without suction. 
^Separation point. 
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TABLE 3 


* 

the FUNCTION G( k , k i) foe the integration 

OF THE DIFFERENTIAL EQUATION 
OF THE MOMENTUM THICKNESS 



> 
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TABLE 3 " Concluded 

THE FUNCTION G( k , tq.) FOR THE INTEGEEIATION 
OF THE DIFFERENTIAL EQUATION OF THE 
MOMENTUM THICKNESS - Concluded 


K-^ e 0 » 3 


X 1 

' 

K 

G('K,«ci) 

0.590 

0.1292 

- 0.542 

.600 

.100 

-.400 

• 6L1 

.0722 

-.265 

.622 

.0463 

-.136 

•637 

.0222 

.0192 

-.016 

0 

.652 

0 

•099 

.688 

-.0380 

.294 

.738 

-.0660 

.438 

.808 

-.0830 

.520 

.892 

-.0910 

•555 

1.172 

-.0721 

.358 


ic i * 0.4 


- 1.099 


X 1 

K 

gOs*].) 

0.770 

0.1068 

- 0.460 

.782 

•0777 

-.324 

.800 

• 0500 

-.200 

.815 

.0202 

.008 

-.080 

0 

.835 

0 

.035 

•875 

-.0415 

.232 

•935 

-.0730 

-.0721 

•375 

.158 


= 0.5 


X 1 


0-935 0.113 
•950 .0825 

.965 .0538 

.982 .0220 

1.0 0 

1.045 -.0460 
1.115 -.0810 

-.0721 


ic ■« *> 0* 6 


G( Kj ) 


MI*1 


-.042 


G( 


gCic,^) 


- 0.125 

-.082 

-.005 

.089 

.262 

.402 






























* 


f 


4 


4 




TABLE 4 


INITIAL VALDES OF TEE PLANE STAGNATION POINT FLOW WITH SUCTION 


H 

X 

K 

g 

«o 

K lo 

C o 

8* 

fl 



C 7 ~ 

hfi 

T 0 

pU | fujiy 

0 

0-355 

-0-645 

0.447 

0.0709 

0 

0 

2-37 

0.266 

0.631 

0.595 

1.163 

.2 

.262 

-.600 

.452 

*0536 

.0904 

• 390 

2.33 

.232 

.540 

• 513 

1-39 

.4 

■177 

"■530 

• 459 

.0372 

_ m 

.104 

•953 

2.28 

•193 

.440 

1 0 rt 

• mzu 

n frO 

4-. (O 

• 5 

*139 

-•473 

.463 

.0299 

.231 

1-345 

2.25 

•173 

• 389 

•374 

2.07 

.6 

.100 

-.420 

.468 

.0219 

.281 

1.900 

2.22 

• 148 

.329 

.316 

2-53 

■7 

*o69 

-•346 

.474 

•0155 

.332 

2.67 

2.18 

•125 

.271 

.264 

3-16 

.8 

.042 

-.256 

.481 

. .0098 

.385 

3*89 

2.12 

•099 

.210 

.206 

4.26 

1.0 

0 

0 

• 500 

0 

• 500 

00 

2 

0 

0 

0 

00 


\J1 


H 
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TABLE 5 


THE BOUNDARY -LAYER PARAMETERS AT THE PLANE PLATE 
IN LONGITUDINAL PLOW WITH HOMOGENEOUS SUCTION; 
LENGTH OF GROWING BOUNDARY LAYER 


1 

X 1 

-K 

&•* 

Sl 

5* 

$ 

To 8 * 

T 0 6* 

T o 

ft 

V 


-PV 0 U 0 

0 

O.OOOITI 
. 000662 
. 00154 
.00291 

0 

.02 

.04 

.06 

.08 

1 

•989 
•979 
• 968 
•956 

1.090 

1.089 

1.088 

1.087 

1.086 

2.66 

2.65 

2.64 

2.63 

2.62 

0 

.0218 

•0435 

.0652 

.0869 

0.572 

•576 

.581 

.586 

•591 

00 

26. 43 
13-34 
8.98 
6.80 

0 

.0131 

.0257 

.0392 

.0540 

• 00*1-56 
.01139 
. 02037 

• 0341 

• 0517 

.10 

.15 

.20 

•25 

.30 

• 945 
.916 
.884 
.851 
.817 

1.085 

1.083 

1.080 

1.077 

1.074 

2.61 

2.58 

2-55 

2.53 

2.56 

.1085 

.1624 

■2159 

.2692 

.3221 

•597 
• 610 
. 625 
.640 

. 656 

5-50 

3-76 

2.89 

2.38 

2.04 

.0676 

.1068 

.1426 

.1845 

.227 

• 0783 
.1124 

• 1551 
.2127 
.2879 

•35 

.40 

•45 

• 50 

• 55 

.780 
.741 
.700 
.656 
. 6lo 

1.070 

1.067 

1.063 

1.059 

1.055 

2.47 

2.44 

2.41 

2.37 

2.34 

• 3746 
.4267 
.4784 
.5296 
.5802 

• 673 
.690 

• 707 
.728 

• 749 

1.80 

1.617 

1.481 

1-375 

I.290 

.280 

•335 
• 394 
.461 
•536 

.3883 

• 5209 

• 7091 

• 9756 
1.373 

.60 
•65 
• 70 
•75 
.80 

.560 

• 507 
.450 
.389 

• 323 

1.050 

1.046 

1.041 

1.035 

1.029 

2.31 

2.27 

2.24 

2.20 

2.16 

.6303 

• 6797 

• 7284 
•7763 
.8233 

• 770 
•793 
.817 
.843 
.871 

1.222 

1.167 

1.122 

1.086 

1.058 

.624 

.722 

.842 

.988 

1.172 

2.007 

3.163 

5-840 

10.556 

14.733 

00 

• 85 

• 90 
•95 
.98 
•99 

1.00 

.252 
.175 
.091 
• 037 
.019 

0 

1.023 

1.016 

1.008 

1.003 

1.002 

1 

2.12 

2.08 

2.04 

2.02 

2.01 

2 

. 8693 

• 9142 
•9578 
■9833 
•9917 

1 

• 900 
•931 

• 964 
.985 

• 993 

1 

1-035 

1.018 

1.007 

1.002 

1.001 

1 

1.416 

1.780 

2.415 

3.250 

3-835 

00 
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TABLE 6 


DRAG LAW OF THE PLANE PLATE IN LONGITUDINAL FLOW 
WITH HOMOGENEOUS SUCTION 


X 1 

S = f(Al) 

F(Ai) 

gUi) 

0 

0 

0 

CO 

.01 

.0000406 

.004665 

115.02 

.02 

. 0001706 

.008883 

52.08 

•03 

.0003682 

.012732 

3^.56 

.014- 

.0006612 

.016623 

25.14 

.06 

.001535 

.02646 

17.29 

.08 

.002909 

.03608 

12.40 

• 10 

. 00456 

.04719 

10-35 

• 15 

.01139 

.07427 

6-528 

.20 

.02037 

.10506 

5.159 

.25 

•03405 

.13949 

4.096 

.30 

.05172 

.18066 

3-493 

•35 

.07833 

.23006 

2.937 

.40 

.1124 

.2872 

2.556 

•45 

•1551 

• 3539 

2.281 

• 50 

.2127 

.4357 

2.048 

• 55 

.2879 

•5357 

1.861 

• 6o 

.3883 

.6614 

1.704 

•65 

•5209 

.8199 

1-574 

• 70 

• 7091 

1.0347 

1.459 

• 75 

•9756 

1-3284 

1.362 

.80 

1.3731 

1.7538 

1.277 

.85 

2.0075 

2.4165 

1.204 

.90 

3.1630 

3.6018 

1.139 

• 95 

5-8403 

6.3098 

1.080 

•98 

10.556 

11.043 

1.046 

• 99 

14.73 

15-23 

1.034 

• 992 

16.15 

16. 64 

1.031 

• 994 

18.01 

18.51 

I.0275 

.996 

20.69 

21.19 

1.0241 

•997 

22.62 

23.12 

1.0220 

• 998 

25.37 

25.87 

1-0196 

♦ 999 

30.11 

30.61 

1.0166 

• 9995 

34.90 

35-40 

1. 0143 

1 

00 

00 

1 


Cf 

■*■00 


= G(Vl) 


c f« = 



Uo 
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TABLE 7 


PARAMETERS FOR THE VELOCITY - DISTRIBUTION OVER THE LENGTH OF GROWING 
BOUNDARY LAYER FOR THE PLANE PLATE IN LONGITUDINAL FLOW 
WITH HOMOGENEOUS SUCTION (TO FIG. 12) 


ft 

X 1 

-K 

f F 

X 1 

• -K 

0 

0 

1 

1.0 

0.754 

0.384 

.1 

.143 

•919 

1.4 

.848 

.255 

.2 

.267 

.840 

1.8 

• 902 

.171 

.4 

• 453 

.696 

3.0 

•973 

.053 

.6 

• 590 

•573 

00 

1 

0 

.8 

.685 

.467 







BACA UM No. 1216 


TABLE 8 


RESULTS OF THE BOUNDARY LAYER CALCULATIONS 
FOR THE CIRCULAR CYLINDER WITH SUCTION 
(a) Cq = 0 
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TABLE 8 - Continued 


RESULTS OF THE BOUNDARY -LAYER CALCULATIONS - Continued 


(B) C Q = 0.5 


_o 1 

<p 

s 

R 


m 

K 

*1 

X 

X 1 

-K 

■8 

0 

0 

0.1573 

0.3681 


0.1112 

0.240 

0.250 

0.580 

2.34 

4 

.06(98 

•1575 

. 3686 

• 0495 

• 1114 

.240 

.250 

.580 


8 . 

•1396 

• 1587 

• 3714 

.0494 

.1122 

.240 

.250 

.580 ! 


12 

.2094 

.1598 

• 3739 

.0493 

.1130 

.240 

.250 

.581 


16 

• 2793 

.1605 

•3756 

.0492 

•1137 

.240 

.250 

.582 


20 

• 349 

.1619 

• 3788 

. 0490 

• 1145 

.240 

.250 

• 584 

2.342 

25 

• 436 

.1634 

.3824 

.0487 

•1155 

.238 

.252 

.585 


30 

.524 

.1658 

.3880 

.0482 

.1172 

.234 

• 255 

•587 


35 

• 611 

•1695 

• 3966 

.0471 

• 1199 

.225 

.2 60 

• 589 


4o 

• 698 

• 1726 

.404 

.0456 

.1220 

.220 

.267 

• 592 

2.343 

45 

• 785 

•1766 

.413 

.0441 

.1249 

.215 

.275 

• 595 

50 

•873 

.1821 

.426 

.0430 

.1288 

.208 

.290 

.601 


55 

•959 

.1871 

.440 

.0403 

•1323 

•195 

.297 

.608 

1 

6o 

•1.047 

• 1937 

• 455 

.0375 

•1370 

.185 

•305 

.615 

i 

65 

1-134 

.2012 

• 473 

• 0331 

.1422 

.164 

.315 

.622 

2.35 

70 

1.222 

.2097 

• 495 

.0287 

.1481 

.141 

.329 

.630 

75 

1.309 

.2181 

• 515 

.0230 

.1542 

•115 

.350 

.640 


8o 

1-396 

.2289 

• 543 

.0170 

• 16L9 

.089 

.360 

. 665 

2.37 

85 

1.484 

.2470 

• 593 

.0109 

.1747 

• 057 

.375 

• 695 


90 

1-571 

.2683 

.645 

0 

• 1897 

0 

.410 

•735 

2.40 

95 

1.658 

.2881 

• 709 

-.0144 

.2037 

-.075 

.462 

.781 


100 

1.746 

.3123 

: .835 

-.0338 

.2208 

-.185 

• 525 

.884 

2.67 

105 

1.833 

• 3421 

• 978 

-.0606 

.2419 

-.415 

.610 

1.06 


S 106. 4 

1.854 

.3522 

1.004 

-.0682 

.2490 

-.452 

.625 

1-18 

2.84 
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TABLE 8 

* 


RESULTS OF THE BOUNDARY Li 

(c) 


9° 

s 

5 

« ./By* 

if v 

5* /UoR 
R 1 V 

n 

0 


0. 3071 


.0698 


.3071 

■O 

• 1396 

• 1351 

.3074 


.2094 

.1360 

.3094 


• 3793 

.13 fir 

•3109 

HS 

•3^9 

.1373 

.3124 

mEm 

.436 

.1383 

• 3146 

30 

.524 

.1400 

.3178 

35 

.611 

.1419 

.3221 

40 

.698 

.1442 

.3273 

45 

.785 

.1466 

.3328 

50 

.873 

.1491 

.3385 

55 

• 959 

.1532 

.3478 

60 

1.047 

• 1595 

•3613 

65 

1.134 

.1665 

.3762 

70 

1.222 

.1729 

.380 

75 

1.309 

.1795 

•394 

8o 

1-396 

.1869 

.410 

85 

1.484 

.1942 

•435 

90 

1-571 

.2017 

.450 

95 

1.658 

.2090 

.476 

100 

1.746 

.2207 

•507 

105 

1.833 

.2291 

.539 

no 

1.920 

.2500 

•598 

115 

2.008 

• 2775 

• 700 

S 115.5 

2.016 

.2814 

•715 


t 
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- Continued 


MLR CALCULATIONS - Continued 
C 0 « 1 


K 

*1 

X 


-K 

0.0365 

0.1909 

0.170 

0.415 

0.525 

•0364 

.1909 

.170 

.415 

.525 

.0362 

.1911 

.170 

.415 

.525 

.0362 

.1923 

.170 

.415 

.525 

.036L 

.1933 

.169 

.417 

.525 

.0354 

.1942 

.166 

.420 

•525 

.0347 

.1956 

.162 

.425 

.524 

.0339 

.1 966 

.156 

.430 

•523 

.0330 

.2007 

• 153 

.440 

.522 

.0319 

.2039 

.150 

.445 

• 521 

.0304 

.2073 

.142 

• 455 

.520 

.0289 

.2106 

•135 

.470 

•519 

.02 69 

.2167 

.127 

.485 

•517 

.0254 

.2256 

.122 

• 495 

.515 

.0234 

• 2355 

.107 

• 512 

• 509 

.0205 

.2445 

.098 

•531 

• 500 

.0167 

.2538 

• 077 

• 554 

• 500 

.0121 

.2643 

•055 

.580 

.508 

.0066 

.2746 

.031 

.600 

•515 

0 

.2852 

0 

.630 

.520 

-.0076 

.2956 

-.040 

.668 

•535 

-.0169 

.3121 

-.088 

• 710 

.560 

-.0272 

.3240 

-.140 

• 775 

•592 

-.0438 

•3^36 

-.225 

.825 

.680 

-.0651 

• 3642 

-•357 

.850 

.835 

-.0682 

• 3697 

-•375 

.870 

.870 
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TABLE 8 


RESULTS OF THE BOUNDARY L 

(a) 


<P° 

e 

R 

m 

Em 


0 

0.1030 

0.2281 


• 0698 

.1030 

.2281 


•1396 

.1031 

.2284 


.2094 

.1034 

.2290 

16 

•2793 

.1039 

.2301 

20 

•349 

.1044 

.2312 

25 

.436 

.1050 

.232 6 

30 

.524 

•1055 

.2340 

35 

.6L1 

.1068 

.2360 

4o 

• 590 

.1078 

.2382 

45 

•785 

.1091 

.2400 

50 

.873 

.1106 

.2437 

55 

• 959 

.1123 

.2472 

6b 

1.047 

•H53 

.2524 

65 

1.134 

.1187 

.2580 

70 

1.222 

.1217 

.2647 

75 

1.309 

.1252 

.2712 

8o 

1.396 

.1289 

.2795 

85 

1.484 

.1334 

.2881 

90 

1.571 

•1393 

.3002 

95 

I.658 

.1460 

.3132 

100 

1.746 

• 1543 

.3290 

105 

1.833 

.1631 

• 3464 

110 

1.920 

.1721 

.3660 

115 

2.008 

• i860 

.3891 

120 

2.095 

.2015 

.415 

125 

2.183 

.2205 

.452 

S 127.5 

2.227 

.2371 

.472 


NACA TM No. 1216 


- Continued 


AYER CALCULATIONS - Continued 

Cq — 2 


K 

D 

a 

H 

-K 

0.0212 

0.2913 

0.090. 

0.620 

0.420 

.0211 

.2913 

.090 

.620 

.420 

.0210 

.2918 

.090 

.620 

.420 

.0208 

.2925 

.087 

.621 

.419 

.0205 

.2939 

.083 

.623 

.418 

.0198 

.2953 

.079 

.625 

.416 

.0194 

.2970 

.075 

.628 

.414 

.0190 

.2984 

.073 

.632 

.412 

.0185 

• 3021 

.071 

.638 

.410 

.0178j 

.3049 

.069 

.645 

.405 

.0168 

.3086 

.068 

.654 

• 399 

.0157 

.3128 

.066 

.665 

• 394 

.0145 

• 3176 

.064 

• 677 

.388 

.0133 

• 3261 

.060 

.690 

.380 

.0119 

•3357 

.052 

.710 

.370 

.0101 

.3442 

.045 

• 735 

.360 

.0081 

•3541 

• 035 

• 755 

.350 

.0058 

• 3646 

.025 

.775 

• 340 

.0031 

•3773 

.012 

• 795 

.325 

0 

• 3940 

0 

.815 

• 310 

-.0040 

.4129 

-.020 

.845 

.295 

-.0083 

.4364 

-.040 

.890 

.278 

-.0138 

• 4613 

-.068 

.925 

•255 

-.0202 

.4868 

-.100 

.980 

.225 

-.0306 

.526 

-.140 

1.055 

.190 

-.o4o6 

• 570 

-.185 

1.190 

.120 

-.0557 

.635 

-.220 

1.260 

.038 

-.0682 

.671 

-.305 

1.305 

-.020 












NACA TM No. 1216 


59 


TABLE 8 - Concluded 


RESULTS OF THE BOUNDARY -LAYER CALCULATIONS - Concluded 


(e) Telocity Distribution 


<p° 

s 

R 

Uk 

cLTJ 

R m 
ds 

0 

0 

0 

2 

4 

.0698 

• 140 

1-995 

8 

• 1396 

.278 

1.981 

12 

.2094 

.416 

1-956 

16 

• 2793 

•551 

1.932 

20 

• 349 

.684 

1.879 

25 

.436 

.845 

1.813 

30 

• 524 

1 

1.732 

35 

.611 

1.147 

1.638 

40 

.698 

1.286 

1.532 

45 

• 785 

1.414 

1.414 

50 

• 873 

1.532 

1.286 

55 

• 959 

1.638 

1.147 

6o 

1.047 

1.732 

1 

65 

1.134 

1.813 

.845 

TO 

1.222 

1.879 

.684 

75 

1.309 

1.932 

.518 

8o 

1.396 

1.970 

• 347 

85 

1.484 

1.992 

.174 

90 

1-571 

2 

0 


f 
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TABLE 9 

RESULTS OF THE BOUNDARY -LAYER CALCULATION FOE THE SYMMETRICAL 
JOUKOWSKY PROFILE J 015 FOR c a = 0 
(a) f x ( o) - 0; C Q = o 


<p 

s 

t 


WM 

K 

K 1 

X 

X 1 

-K 

180 

0 

0.0370 

0.0878 

0.0708 

0 

0.355 


D 

0.650 

177-5 

-0493 

.0377 

.0895 

• 0697 



.349 



•653 

175 

. 00986 

.0405 

.0964 

.0674 



*335 



.661 

172-5 

-0148 

.0438 

.1048 

.0630 



• 319 



.682 

170 

-0197 

.0485 

.1167 

.0581 



.298 



• 703 

165 

-0308 

.0605 

.1480 

.0464 



.247 



•753 

160 

.0444 

-0755 

.1885 

• 0351 



.197 



.802 

150 

-0764 

.1109 

.284 

• 0193 



• 108 



.892 

140 

.1233 

-1572 

.411 

• 0053 



.030 



-964 

136 

.1445 

.1772 

.472 

0 



0 



1.000 

130 

-1775 

.2062 

.559 

-.0089 

■ 

■ 

-.052 



1.062 

120 

.2416 

.2608 

.742 

-.0260 

I 

■ 

-.172 



1.182 

no 

.3132 

.3217 

.960 

-.0484 

1 

■ 

-.362 


J 

1.367 

s 101-9 

• 377 

.3748 

1.310 

-.0680 

B 

■ 

-.628 


f 

1.637 
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TABLE 9 " Continued 

■¥ 


RESULTS OF THE BOUimAItf -LAYER CALCULATION - Continued 


(1)) f^o) = 0.5; Co - 0.0695 


<fp 

B 

t 

?/? 

5»^o^ 

K 

K 1 

X 

X 1 

-K 

l8o 

0 

0.0359 

0.0846 

0.0667 

0.0180 

0.335 

0.044 

0.640 

177*5 

-.00493 

.0366 

.0867 

.0658 

.OI83 

*331 

.045 

.645 

175 

.00986 

*0391 

.0927 

.0634 

.-0196 

.322 

• 046 

.654 

172*5 

.0148 

.0430 


.0603 

.0215 

•307 

.047 

.674 

170 

*0197 

.0474 

.n4o 

*0556 

.0237 

.283 

.051 

.690 

165 

.0308 

.0580 

.1416 

.0430 

.0299 

.227 

.060 

.7^4 

l6o 

.0444 

*0737 

.1830 

.0334 

.0368 

.180 

.074 

*773 

150 

. .0764 

.1072 

.268 

.0169 

.0540 

.100 

.110 

•833 

l4o 

•1233 

.1446 

*367 

.0045 

.0720 

.024 

.154 

.867 

136 

.1445 

*16L3 

.415 

0 

• 0790 

0 

*173 

.880 

130 

*1775 

.1850 

.483 

-.0071 

.0913 

-.035 

.203 

.906 

120 

.2416 

.2274 

.620 

-.0200 

-1137 

-.112 

• 255 

*976 

110 

.3132 

*2793 

• 764 

-*0356 

*1397 

-.225 

.316 

1*075 

100 

*391 

*3339 

.961 

-.0530 

.1669 

-.360 

.384 

1.200 

S 92.lt- 

*4575 

.'3742 

1.125 

-.0682 

.1871 

-.505 

.440 

1.310 
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TABLE 9 “ Continued 


RESULTS -OF THE BOUNDARY -LAYER CALCULATION - Continued 
(c) f^o) = l.o; C Q = 0.139 


<p° 

s 

l 

tS /u c t 
ij “v“ 

♦1+r 

(Ol 

K 

*1 

X 

X 1 

-K 

180 

0 

0.0348 

0.0819 

0.0626 

0.0348 

0-314 

0.077 

1 

0.631 

177.5 

. 00493 

.0358 

.0842 

.0626 

.0358 

•313 

•079 

.631 

175 

.00986 

.0382 

.0900 

.0605 

.0388 

.308 

.085 

.640 

172.5 

.0148 

• o4i6 

.0985 

.0564 

.0403 

.293 

•095 

•653 

170 

• 0197 

.0460 

.1095 

.0516 

.0450 

.265 

.104 

.667 

165 

.0308 

.0560 

•1345 

■ .0376 

.0560 

.205 

.128 

.703 

160 

.0444 

• 0686 

.1678 

.0292 

.0678 

.154 

.156 

.748 

150 

.0764 

.0985 

.243 

.0152 

• 0985 

.080 

.223 

.776 

140 

.1233 

.1342 

•332 

.0038 

.1340 

.020 

.300 

.788 

136 

.1445 

.1484 

.368 

0 

.1478 

0 

•335 

.795 

130 

•1775 

.1682 

.421 

-.0059 

.1682 

-.032 

.380 

.800 

120 

.2416 

.2030 

.505 

-.0162 

.2025 

-.084 

.465 

■ 806 

110 

.3132 

.2385 

• 595 

-.0270 

.2400 

-.142 

• 556 

.809 

100 

• 391 

.2790 

.696 

-.0380 

.2800 

-.200 

.650 

.812 

90 

• 475 

.3200 

.800 

-.0499 

.3203 

-.267 

.750 

.818 

80 

• 561 

.362 

• 916 

-.0636 

• 3622 

-.343 

.858 

.830 

S 76.4 

-593 

• 376 

• 952 

-.0682 

•3766 

-.380 

.892 

.860 
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TABLE 9 " Continued 


RESULTS OF TEE BOUirDAEr -LAYER CALCULATION - Continued 
(d) fj^O) = 1-55 C 0 - 0.2085 




<p° 

s 

t 

mm 

B 

K 

K 1 

X 

X 1 

-K 

180 

0 


0.0792 

0.0584 

0/0504 

0.295 

0.110 

0.635 

177.5 

.00493 


.0824 

.0592 

.0526 

• 293 

.116 

.635 

175 

.00986 

.0372 

.0878 

.0580 

• 0557 

.281 

.125 

.639 

172.5 

.0148 

.0404 

.0958 

■ 0531 

.0606 

.265 

.136 

.648 

170 

.0197 

.0442 

.1052 

.0482 

.0662 

.202 

.150 

.658 

165 

.0308 

.0544 

.1301 

.0376 

.0816 

.192 

.181 

• 690 

160 

.0444 

.0675 

.1596 

.0280 

.0980 

.145 

.220 

•705 

150 

.0764 

.0927 

.224 

•0135 

• 1391 

.070 

.310 

.722 

140 

•1233 

.1249 

.297 

.0027 

.1852 

.016 

.415 

.718 

136 

.1445 

.1368 

.325 

0 

.205 

0 

.454 

.702 

130 

.1775 

.1543 

.365 

-.0046 

.230 

-.014 

•512 

• 669 

120 

.2416 

.1825 

.430 

-.0127 

.272 

-.060 

.605 

.630 

no 

.3132 

.2107 

.490 

-.0202 

•315 

-.097 

.698 

• 596 

100 

•391 

.2388 

• 546 

-.0276 

• 359 

-.132 

.785 

.566 

90 

.475 

.2683 

.610 

-.0356 

.401 

-.165 

.870 

.529 

80 

.561 

.2981 

. 668 

-.0432 

.443 

-.198 

• 950 

.486 

70 

.645 

.3255 

.726 

-.0504 

.487 

-.230 

1.034 

.430 

60 

.734 

•3555 

• 770 

-.0578 

•531 

-.260 

1.114 

•350 

50 

.807 

.3789 

.808 

-.0637 

.568 

-.286 

1.190 

.250 

S 40 

.885 

.4027 

.835 

-.0682 

.604 

“•305 

1.265 

.150 
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TABLE 9 ” Continued 


RESULTS OF THE B OUNDAKY -LAYER CALCULATION - Continued 
(e) fi(O) =3-0} C Q = 0.417 


<p° 

£ 

t 


gg 

K 

*1 

X 

h . 

-K 

l8o 

0 

0.0307 

0.0720 

0.0504 

0/0921 

0.242 

0.210 

0.602 

177-5 

.00493 

.0323 

.0755 

.0511 

.0940 

.241 

.213 

.603 

175 

. 00986 

.0342 

.0795 

.0494 

.1026 

.229 

.225 

.604 

172-5 

.0148 

.0365 

.0852 

.0433 

.1095 

.212 

.240 

.609 

170 

• 0197 

.0398 

.0900 

• 0391 

.1194 

• 195 

.2 60 

.612 

165 

.0308 

.0486 

.1150 

.0300 

.1398 

• 157 

.310 

.620 

160 

.0444 

.0587 

.136 

.0216 

•1705 

.107 

■376 

. 6lo 

150 

.0764 

.0794 

.184 

.0099 

.238 

.051 

.530 

• 570 

140 

■1233 

.1015 

.231 

.0022 

• 309 

.oil 

.666 

.502 

136 

.1445 

.1118 

.248 

0 

•335 

0 

• 714 

.471 

130 

•1775 

.1240 

.270 

-.0032 

• 370 

-.016 

.781 

.422 

120 

.2416 

.1435 

.303 

-.0080 

.425 

-.033 

. 880 

.306 

110 

.3132 

• 1594 

•333 

-.0119 

.476 

-.049 

• 970 

•153 

100 

• 391- 

• 1745 

• 358 

-.0148 

.522 

-.060 

1.045 

.022 

90 

• 475 

.1881 

• 374 

-.0172 

• 564 

-.068 

1.110 

-.090 

80 

• 5 6l 

.2006 

• 396 

-.0197 

.600 

-.072 

1.172 

-.180 

70 

.645 

.2109 

.410 

-.0212 

•633 

-.075 

1.230 

-.230 

60 

• 734 

.2218 

.424 

-.0226 

.663 

-.078 

1.280 

-.270 

50 

.807 

.2302 

• 439 

-.0234 

.691 

-.079 

1.326 

-.300 

40 

.885 

.2383 

.452 

-.0238 

• 715 

-.080 

1-336 

-.330 

30 

• 945 

.2440 

.464 

-.0248 

..732 

-.080 

1-340 

-•350 

20 

•990 

.2486 

.470 

-.0256 

• 746 

-.080 

1.344 

-.360 

10 

1.018 

.2512 

.476 

-.0259 

•754 

-.080 

1-348 

-.360 

0 

no sepai 

1-028 
•at ion 

.2522 

• 479 

-.0261 

• 757 

-.080 

1-350 

-.360 
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* TABLE 9 " Concluded 

RESULTS OF THE BOUNDARY -LATER CALCULATION - Continued 
(f) "Velocity Distribution 


<p° 

s 

t 

U 

U 0 

t dU 
U 0 ds 

180 

0 

0 

51.7 

ITT- 5 

.00493 

.231 

49.1 

175 

.00986 

.445 

42.2 

172.5 

.0148 

.632 

32.6 

170 

.0197 

.782 

24.7 

165 

.0308 

.993 

12.69 

160 

.0444 

1.111 

6.16 

150 

.0764 

1-233 

1.572 

140 

-1233 

1.270 

.215 

136 

.1445 

1.271 

0 

130 

.1775 

1.267 

-.209 

120 

.2416 

1.245 

-.387 

110 

.313 

1.212 

-.468 

100 

• 391 

1.178 

-.485 

90 

• 475 

1.136 

-.487 

80 

.561 

1-095 

-.487 

70 

.645 

1-053 

-.475 

60 

.734 

1.014 

-.458 

50 

.807 

• 977 

-.442 

40 

.885 

■ 945 

-.420 

30 

•945 

• 917 

-.416 

20 

• 990 

.899 

-.414 

10 

1.018 

.889 

-.410 

0 

1.028 



CO 

• 

-.410 


« 
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TABLE 10 

VELOCITY DISTREBUTION OVER TEE REGION CF GROWING BOUNDARY LAYER 
FOR THE ELATE WITH HOMOGENEOUS SUCTION 


ft 

0.1 

(T 

= 0.2 

it - 

» 0.4 

VO 

6 

0 

f T 

«= 0* 8 

-▼oy 

u 

'W 

U 


U 


U 


V 

V 

% 

V 

Uo 

V 

*0 

V 

u 0 

V 

u 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.0286 

.1108 

•0534 

.1168 

.0910 

•1279 

.1180 

-1373 

•137 

•1455 

• 0572 

;2l76 

.1068 

.2271 

.1820 

.2448 

.236 

.2598 

.274 

.2727 

.0858 

.3205 

.1602 

•3317 

.273 

.9522 

.354 

• 3697 

.411 

• 3848 

.11 44 

-4173 

.2136 

.4287 

."*64 

•4497 

.472 

.4676 

• 548 

.4830 

• 143 

•5107 

.267 

.5211 

.455 

• 5302 

.590 

.5564 

.685 

■5704 

.1716 

-59 #7 

.3204 

.6054 

• 546 

.6215 

.708 

• 6351 

.822 

.6469 

.2002 

.6760 

.3738 

.6826 

.637 

.6948 

.826 

•7051 

• 959 

• 7141 

.2288 

•7475 

.4272 

• 7518 
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The parameter K and (See table 70 
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Figure 3. - The form parameter K of the velocity profile as a function of 

A , Ap according to equation (14). 
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Figure 8.- Comparison of the velocity distributions according to the 
approximate with the exact calculation. 

(a) Plane plate in longitudinal flow, exact calculation according to Blasius, 

approximate calculation according to equations (53) and (54). 

(b) Exact calculation according to Hiemenz, approximate calculation 

according to equations (58) and (58a). 
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Figure 9.- The extent of growing boundary layer for the plane plate with 

-v 0 B* 5 * t B* . r — - 

homogeneous suction: — - — * — * — — against \J s . 


J 
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Figure 10.- Plane plate with homogeneous suction; region of growing 
boundary layer; velocity distribution. 
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Figure 11.- Plane plate with homogeneous suction; region of growing boundary 

T 0 j— 

layer; the local friction coefficient against v S . 

O oo 






Figure 13.- Plane plate with homogeneous suction, region ol growing boundary layer: The coefficient of the 
total drag c f ag ains t U 0 l/v for various values of suction coefficient -v 0/ /U Q . 
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Figure 14.- Plane stagnation point flow: comparison of the boundary^layer 
thickness and momentum thickness of approximate and exact calculation 
for various suction quantities, exact calculation according to reference 9. 
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Figure 15.- The boundary layer on the circular cylinder with homogeneous 

suction for various suction quantities C n = \ — 8 — . Form 

2 U ° V 
parameter k = With increasing suction quantity the 

separation point shifts rearward. 
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Figure 16.- The boundary layer on the circular cylinder with homogeneous suction for the mass 
coefficient C Q = 0.5. Comparison of the approximate calculation with an exact solution 

( r o-foi'onr’ o 17^ 
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Figure 17.- The boundary layer on a symmetrical Joukowsky profile J 015 
with c a = 0 for homogeneous suction with various suction quantities 

Cr, = according to equation (36): C_ = — x f,(o) 

' n/kT 2 & 

(K^ = 51.7). R = ~ . With increasing suction quantity the 

separation point shifts rearward. 








' ^ oe a* f o~e o.b y^r 10 

V % 

Figure 19, - Region of growing boundary layer for the plane plate with homogeneous suction; comparison of 

approximate and' exact calculation for the wall shearing stress. 
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